Chin. J. Astron. Astrophys. Vol.0 (200x) No.O, 000-000 Chinese Journal of 

[http://www.chjaa.org ) Astronomy and 

Astrophysics 



Effects of Magnetic Fields on Neutrino- dominated 
Accretion Model for Gamma-ray Bursts 

Yi Xie*, Chang- Yin Huang and Wei-Hua Lei 

Department of Physics, Huazhong University of Science and Technology, Wuhan 430074, 
China 



Abstract Many models of gamma-ray bursts suggest a common central engine: a 
black hole of several solar masses accreting matter from a disk at an accretion rate 
from 0.01 to 10 M Q s _1 . The inner region of the disk is cooled by neutrino emission 
and large amounts of its binding energy were liberated, which could trigger the 
fireball. We improve the neutrino-dominated accreting flows by considering the 
effects of the magnetic fields, and find that more than half of the liberating energy 
can be extracted directly by the large-scale magnetic fields on the disk. And it turns 
out that the temperature of the disk is a bit lower than the neutrino-dominated 
accreting flows without magnetic field. Therefore, The outflows are magnetically- 
dominated rather than neutrino dominated. In our model, neutrino mechanism can 
fuel some GRBs (not the brightest ones), but cannot fuel X-ray flares. However, 
the magnetic processes (both BZ and electromagnetic luminosity from a disk) are 
viable mechanisms for most of GRBs and the following X-ray flares. 
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1 INTRODUCTION 

Gamma-ray bursts (GRBs) are flashes of gamma-rays occurring at cosmological distances, 
being the most powerful explosions since the Big Bang. They are generally divided into two 
classes (Kouveliotou et al.1993): short-duration (Tgo < 2s) hard-spectrum GRBs (SGRBs) 
and long-duration (Tgo > 2s) soft-spectrum GRBs (LGRBs), which have different progenitors. 
LGRBs root in core collapses of massive, rapidly rotating stars (Woosley 1993, Paczynski 1998, 
Hjorth et al 2003, Stanek et al 2003), and supernovae have been observed coincidently in some 
LGRBs (Galama et al 1998; Stanek et al 2003; Hjorth 2003). In contrast to LGRBs, SGRBs 
may arise from coalescence of neutron stars or black hole binary systems due to damping of 
gravitational radiation (e.g. Eichler et al. 1989, Narayan, Paczynski & Piran 1992, Fryer & 
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Woosley 1998), and they are probably associated with elliptical galaxies (Gehrels et al. 2005; 
Bloom et al 2006; Barthelmy et al. 2005; Berger et al 2005). It is believed that the two processes 
give rise to a black hole of several solar masses with a magnetized disk or a torus around it 
(Meszaros & Rees 1997b). And many central engine models of GRBs based on this scenario 
(exception models, for instance, magnetizedrotating neutron stars, see e.g. Usov 1992). 

Some authors have studied the accretion model for GRBs by assuming steady-state accretion 
(e.g. Papham, Woosley & Fryer 1999, hereafter PWF; Narayan, Piran & Kumar 2001, hereafter 
NPK; Di Matteo, Perna & Narayan 2002, hereafter DPN). Their studies show that at the 
extremely high accretion rate (0.01 to 10 M s _1 ) needed to power GRBs, the disk cannot be 
cooled efficiently as the gas photon opacities are very high, and a large fraction of its energy is 
advection dominated. However, inner region of the disk becomes hot and dense enough to cool 
via neutrino emission, and this accretion mode is referred to as neutrino-dominated accretion 
flows (NDAFs). The neutrinos can liberate large amounts of binding energy via the vv — > e + e~ 
processes in regions of low baryon density and then trigger the fireball. 

However, the model with "neutrino-driven outflow" alone cannot be a candidate of some 
GRBs central engine. For instance, numerical simulations by Shibata et al. (2006) suggest that 
the collapse of hypermassive neutron-star triggered by gravitational wave cannot be a candidate 
for the central engine of SGRBs, however, it becomes powerful enough to produce the fireball 
after taking the magnetic braking and MRI into account. On the other hand, researches show 
that the magnetic fields can be magnified up to 10 15 ~ 10 16 G by virtue of magnetorotational 
instability (MRI, Balbus & Hawley 1991) or dynamo processes (Pudritz & Fahlman 1982 and 
references therein) in the inner region of the disk. So, the existence of strong magnetic fields 
should be considered. Both PWF and DPN compared the luminosity of neutrino emission and 
Poynting flux, and indicated that MHD processes are viable mechanisms for powering GRBs, 
but they did not include magnetic fields in their disk conditions. 

In this paper, we intend to improve the NDAF model by considering the effects of magnetic 
braking and magnetic viscosity. The equation of angular momentum of a standard disk is 
replaced by the equation of a magnetized disk in which the viscosity caused by magnetic braking 
and magnetic viscosity only. Meanwhile, we deduce the rotational energy extracted by large-scale 
fields in the disk from the thermal energy produced by viscous dissipation, and the magnetic 
pressure is considered in equation of state else. It turns out that the inner region of the disk 
is magnetically dominated. Magnetized accretion models within the GRB context have been 
also discussed in several other papers in which detailed numerical simulations are performed. 
For example, Proga et al. (2003) studied a MHD collapsar accretion model and suggested that 
MHD effects alone can launch a GRB jet, which is Poynting flux dominated. Mizuno et al. 
(2004a; 2004b) drew similar conclusions using a GR-MHD code. These results agree with our 
conclusion in this paper. 

This paper is organized as follows. In Section 2 we outline the theory of a magnetized 
accretion disk. In Section 3 we introduce the basic assumptions and equations of our model. In 
Section 4 we show the numerical results of our model, and finally, In Section 5 we summarize 
main results of our model and some defects are discussed. 

2 DESCRIPTION OF A MAGNETIZED ACCRETION DISK 

It is widely known that the magnetic fields on the disk can greatly affect angular momentum 
transfer and hence the accretion rate via a variety of modes. In this paper we only consider 
two basic mechanisms: the first one is magnetic viscosity, the weak magnetic fields creates MRI 
in which turbulence dominates the angular momentum transfer (Balbus & Hawley 1991). The 
second is the magnetic braking, the large scale magnetic field extracts rotational energy of 
disk due to the shear force of differential rotation (Blandford 1976; Blandford & Payne 1982; 
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Livio et al 1999). We assume that the accretion process is governed by these two mechanisms 
completely. The main points of the magnetized accretion disk model (in which the magnetic 
braking and magnetic viscosity are considered only for angular momentum transfer) given by 
Lee et al. (2000) are outlined as follow. 

According to Torkelsson et al. (1996) the magnetic viscosity v ma9 is defined as 

B4.Br/ 4tt = -v mag (rdn dl sk/dr) p, (1) 

where B^ and B r are respectively the azimuthal and radial components, Sldi.sk is angular velocity 
and p is the density of the disk matter. The magnetic viscosity v mag can be parameterized as 
(Shakura & Sunyaev 1973; Pringle 1981) 

v mag = ( 2 ) 
1/2 

where c s is the sound velocity of the disk (c s = (Pdisk/p) , m which Pdisk is the disk pressure) 
and H is the half-thickness of the disk. Invoking hydrostatic equilibrium perpendicular to the 
disk plane, we have H = c s /Sldisk and 

v mag = a ™a g Jj^.^ (3) 
1 /2 

For a Keplerian orbit we have Vldi.sk ~ = (GM/r 3 ) , and Eq. (1) can be written as 

B4,B r /AK = \a ma °P disk . (4) 

The accretion rate is generally determined by magnetic braking for H <C r (Lee, Wijers & 
Brown 2000), the angular momentum balance equation can be written as 

M = 2rB4B z /n dlsk - (5) 
The axisymmetric solution (Blandford 1976) is 

B ljj = 2rVL d iskBz/c. (6) 

A roughly steady state will be reached when the grown rate of B^ generated by differential 
rotation from radial field equals to its loss rate by buoyancy, then the magnitude of B^ can be 
estimated as (Katz 1997) 



B4, 

By using Eqs. (4), (6) and (7) we have 



\B r n disk H 



1/2 

(4^) 1/4 - (7) 



Bz = l(^^\'\^f'\ (8) 



2 V GM 
and 

B4 = (nPMsk) 1/2 (9a ma °) 1/3 . (9) 

The vertical component and azimuthal component of field can be estimated by using Eqs. 
(8) and (9), and for a given a ma9 , only depending on the gas pressure. 

Combining Eqs. (6) and (8) with Eq. (5) we have (see also Lee, Wijers & Brown 2000) 

M = 4r 2 B 2 z /c. (10) 
Since M is independent of r in steady-state accretion, we can infer that Bz oc 1/r. 
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3 BASIC ASSUMPTIONS AND EQUATIONS OF MODEL 

The basic physical conditions in disk models for GRBs can be derived by virtue of steady- 
state conditions (PWF, NPK, DPN). Base on these studies, we consider the effects of magnetic 
field in inner regions of disks in the frame of hydrodynamics. The basic equations consist 
of equation of state and the conservation equations of energy and angular momentum in a 
magnetized accretion disk, which are described as follows. 

In equation of state we include the contributions from radiation pressure, gas pressure, 
degeneracy pressure and magnetic pressure, 

11 ^ 4 , pkT , 2nhc ( 3 V /3 ( p \ i/3 , B 2 



P= aT^ + t— + — [-— ^ +— , (11) 

12 m p 3 \8irmp J \PeJ 87r 

where a is the radiation constant, T is the disk temperature, and the factor includes the 
contribution of relativistic electron-positron pairs. In degeneracy term, p e is the mass per elec- 
tron, and it is taken as 2 by assuming equal number of protons and neutrons. For the magnetic 
pressure we only consider poloidal component in calculation, provided that it is not much less 
than the toroidal component. 

The conservation of mass is written by NPK and DPN as follows, 

M = ATrrv r pH w 6irpvH, (12) 

where v r is the radial velocity and v r = 3v/2r. Different from NPK and DPN, we replace the 
Eq. (12) by Eq. (10), which includes the effects of magnetic braking and magnetic viscosity. 

In energy equation, the viscous heating equals neutrino radiative loss plus advective loss 
and the fraction of rotational energy extracted by large-scale magnetic fields 

2>GMM , _ _ N 

87rr 3 = + Q eN ) H + Qadv + Q B . ( 13 ) 

in which q'^ is cooling via pair annihilation and we take it as the approximation of Itoh 
ct al. (1989; 1990): ~ 5 x 10 33 T^er g s cm' 3 s^ 1 (in which X n = X/10 n is used). The 
q~ N represents the cooling via pair capture on nuclei, and can be estimated as q~ N ~ 
9 x 10 33 /9ioTj 3 1 er(7scTO _3 s _1 . And q a & v is the advective cooling rate, we approximate it by 
(see eg., Narayan & Yi 1994; Abramowicz et al.1995) 

, ds ..A/H.rf . ipkTl + Xnuc 



q adv = H V T-~i V -Tl—aT A + -'- , (14) 

dr r \ 3 2 m p 4 J 

in which s is specific entropy, X nuc is the mass fraction of free nucleons, £ oc —dlns/dlnr 
is assumed to be equal to 1 as in DPN. And finally, Qg represents the energy extracted by 
magnetic field (see, Lee, Wijers & Brown 2000) 

dP ma s _ B 2 z r 2 (GM\ _ GMM 
Qb ~ dS ~ nc { r 3 ) ~ Airr* ' ( 5) 

where dS — 2irrdr. Comparing Eq.(15) with Eq.(13), we find that two thirds of energy of 
viscous heating was substituted by the field extracting, therefore, it is magnetically dominated 
in the inner regions of disks. 

Eqs. (8), (10), (11) and (13) contain four independent unknowns P, p, T and E>z as functions 
of r and compose a complete set of equations which can be numerically solved with given M, 
a (for simplicity, we omit the superscript 'mag') and M. In the following calculations we fix 
M = 3M Q (the corresponding Schwarzschild radius R s is 2GM/c 2 = 8.85 x 10 5 cm), a = 0.1 . 
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Fig. 1 Pressure components profiles for three values of the accretion rate: (a) rh = 0.1, 
(b) to = 1, and (c) rh = 10. The gas pressure is shown by the solid line, degeneracy 
pressure by dotted line, radiation pressure by dashed line, and magnetic pressure by 
long-dashed line. 



4 NUMERICAL RESULTS 

4.1 Gas Profiles 

We show the numerical solutions of the full equations in this section, and the software of 
" Mathematica" is used for the numerical algorithm of Newton iteration method. The pressure 
components profiles are shown in Figure 1, and the solutions for three values of the accretion rate 
rh = 0.1, 1, and 10 (m is defined as rh = M/M Q s _1 ) are show in (a), (b) and (c) respectively. 
The gas pressure, degeneracy pressure, radiation pressure, and magnetic pressure are shown by 
the solid line, dotted line, dashed line, and long-dashed line, respectively. From Figure 1, we 
obtain the following results: 

(i) From (a) we can see that, the flows is radiation pressure dominated in the region of 
IRs ~ 10 Rs and may be thermally unstable (see 4.3, stability analysis). It is thermally stable 
in the same region in DPN as it always dominates by gas pressure. 

(ii) From (b) and (c) we can see that, the magnetic pressure component is more important 
at large radii and even overwhelms the gas pressure and degeneracy pressure. So our model is 
valid only in a narrow region because of the restriction of Eq. (12). 

Temperature and density profiles calculated from our model are shown in Figure 2. We 
show our solutions for three values of the accretion rate, m = 0.1, 1, and 10 (long dashed, solid, 
and short-dashed lines, respectively). Comparing with DPN (see DPN, Fig.l), we find that the 
temperature of disk is a bit lower than NDAFs without considering the effects of magnetic 
fields, and the density drops much more rapidly with the radius. 

4.2 BZ Luminosity, Electromagnetic Luminosity from Disk, and Neutrino 
Luminosity 

It is a common assumption that the magnetic fields will rise up to some fraction, which for 
instance in DPN, 10% of its equipartition value B 2 /8tt ~ pc 2 s . For 0.1 < to < 10, the typical 
values of pc 2 are 10 30 ~ 10 32 ergscm~ 3 , implying a field strength of 10 15 — 10 16 G. The BZ jet 
luminosity is then 
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Fig. 2 Temperature and density profiles in (a) and (b) respectively. The profiles are 
shown for three values of the accretion rate: m = 0.1 (long dashed lines), m = 1 (solid 
lines), and to — 10 (short dashed lines). 



L BZ = %^R% * 10*V (4§c)\^ er 9 scm-*, (16) 

in which a is the dimcnsionlcss black hole spin parameter, Bh is the magnetic field on the 
horizon . 

The electromagnetic power output from a disk is equal to the power of the disk magnetic 
braking can be calculated as (Livio et al. 1999; Lee et al. 2000) 

Ld =^ r {—) cRia UJ Lbz - (17) 

Consistent with previous work (Merloni & Fabian 2002; DPN), we take approximately, 

B z ~(H/r)B H . (18) 

It is easy to get the strength of poloidal field in the disk and the field on the black hole 
horizon for a given to by using Eqs. (10) and (18) in our model, without the assumption of 
equipartition value discussed above. And then, the BZ jet luminosity and the electromagnetic 
power from a disk can be calculated from Eqs. (16) and (17). The neutrino luminosity is given 
by L v = J^ m,n 2nq~rdr, in which q~ = {q~ D + q~ N ) H, r min = IRs (for an extreme Kerr black 
hole), and r m i n = IO-R5. We estimate the luminosity due to vv annihilation along z-axis above 
the disk to be the Eq. (21) in DPN. In Figure 3 we show the curves of BZ luminosity Lbz (solid 
lines), electromagnetic luminosity from a disk Ld ( long-dashed line), and neutrino annihilation 
luminosity L vT> (short-dashed line) versus dimensionless accretion rate. From Figure 3, we obtain 
the following results: 

(i) Ld is larger than Lbz for H/r = 0.2 (solution of this model), however, both of them are 
viable mechanisms for central engines of GRBs , and can also fuel the observed X-ray flares in 
which case the accretion rate of m = 0.01 is needed. 

(ii) The L vi? is around 10 51 ergss _1 at to = 1, which is sufficient to power some GRBs, and 
at the accretion rate of m = 0.01, it fails to fuel the X-ray flares. 
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Fig. 3 The solid lines represent the BZ luminosity, and the long-dashed line represent 
the electromagnetic luminosity from a disk and short-dashed line represents neutrino 
annihilation luminosity as a function of accretion rate and value of H/r = 0.2. 



In conclusion, neutrino mechanism can fuel some GRBs (not the brightest ones), but cannot 
fuel X-ray flares. However, the magnetic processes (both BZ and electromagnetic luminosity 
from a disk) are viable mechanisms for most of GRBs and the following X-ray flares (this 
agrees well with discussions of Fan et al. (2005)). 



4.3 Stability 

Both NPK and DPN discussed the stability properties of their solutions. Since our model 
considers the effects of magnetic fields and differs considerably with theirs, it is interesting to 
examine whether the solution is stable. 

The general condition for thermal stability is given by (Piran 1978) 

/dlnQ+\ (<nnQ_\ 

in which Q ± are the integrated (over the height of the disk) heating (+) and cooling (-) rates. 
The cooling rate Q~ = q~ + q a dv + Qb- We show the two curves of Q~ and Q + as a function 
of gas temperature in Figure 4. The radius is fixed at r = 5i?s, while the surface density 
is taken to be £ = I0 16 gcm~ 2 . From Figure 4 we can see that, the flow is unstable while 
the temperature T is lower than 5 x lO 10 ^, because the magnetic fields extract rotational 
energy from disk is independent with temperature. When T > 5 x lO 10 ^, it turns to stable 
because q~ N oc T 6 becomes relatively significant with the temperature increasing. When the 
disk temperature crosses the critical point of the instability curves, the thermal energy would 
be released suddenly in a thermal time scales. It is possibly an explanation for the variability 
time scales of tens of msec in the light curves of the GRBs, and we will give the details in 
another paper. 

Following NPK and DPN we use the condition for viscous stability 
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Fig. 4 Thermal stability analysis. Dashed line shows Q + and solid line shows Q 
vs. T for r = 5Rs curves, and for X = I0 le gcm~ 2 . The solution is unstable while 
T < 5 x 10 W K , and becomes stable when T > 5 x lO 10 ^. 



In our model, we have M oc £ for the case of gas pressure dominated, M oc £ 8 / 7 for degeneracy 
pressure case, M oc £ 2 for radiation pressure case, and we have M oc £? 2 oc P for the magnetic 
dominated and the surface density S oc P 1 / 2 p 1 / 2 , meanwhile, considering the magnetic fields 
decreases as r -1 , we have approximately p oc r~ 3 , and then M oc E. All of these cases are 
clearly viscously stable. 

Finally, we also consider the gravitational instability. The accretion flow will become grav- 
itational unstable if the Toomre parameter Qt is less than 1, for a Keplerian orbit, which is 
given by (Toomre 1964) 

Qt = ^ = ^- (21) 

We have checked that Qt ^ 1 hence the flow is gravitationally stable in the inner region of 
the disk. Nevertheless, at large radii the Toomre parameter could be less than 1 as argued by 
Perna et al. ( 2006 ). Actually that was another model for X-ray flares. 



5 CONCLUTION AND DISCUSSION 

In this paper we modify the NDAFs model as a central engines for GRBs by considering 
the effects of magnetic braking and magnetic viscosity in the frame of Newtonian dynamics. We 
found that two thirds of the liberating energy was extracted directly by large-scale magnetic 
fields on the disk and the temperature of a disk is a bit lower than the NDAFs without magnetic 
fields. Furthermore, the density of the disk drops faster than NDAFs along the radius. Therefore, 
the inner region of the flow is magnetically dominated rather than neutrino dominated. However, 
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the neutrino mechanism can still fuel some GRBs (not the brightest ones), but cannot fuel X- 
ray flares. However, the magnetic processes (both BZ and electromagnetic luminosity from a 
disk) are viable mechanisms for most of GRBs and the following X-ray flares. 

Our model is formulated invoking Newtonian potential, ignoring the effects of general rela- 
tivity which may be important in some aspects (Gu et al. 2006) and neutrino opacity. Specially, 
the main simplification of the analytic approach (this and other analytical works) is the require- 
ment of a steady state, which in reality (e.g. numerical simulations) is not necessarily justified. 
An example is that with magnetic fields, both numerical simulations (Proga & Begelman, 2003) 
and some analytical arguments (Proga & Zhang, 2006) suggest that the accretion flow may not 
always in a steady state. Rather, magnetic fields accumulated near the black hole can form a 
magnetic barrier that temporarily blocks the accretion flow. This makes some dormant epochs 
at the central engine. The breaking of the barrier would lead to restarting the central engine, 
which is required to explain the recent Swift observations of X-ray flares (for a review of Swift 
results and in particular X-ray flares and their interpretations, see Zhang, 2007). In this pa- 
per, the poloidal component of magnetic fields is Bz oc 1/r, which implies that the magnetic 
pressure drops much slower than the other components and the calculations indicate that the 
magnetic fields pressure could be dominant at larger radii. In fact, such over-pressure mag- 
netic fields are the agent to form the magnetic barrier as reported by Proga & Zhang, which is 
needed to interpret the observed X-ray flares. The unsteady state accretion model and the case 
of over-pressure magnetic fields will be studied in our future work. 
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